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Abstract 

It has been known for some time that 2-loop renormalization group (RG) equations of a dimen- 
sionless parameter can be solved in a closed form in terms of the Lambert W function. We apply 
the method to a generic theory with a Gaussian fixed point to construct RG invariant physical 
parameters such as a coupling constant and a physical squared mass. As a further application, we 
speculate a possible exact effective potential for the O(N) linear sigma model in four dimensions. 
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I. INTRODUCTION 



The purpose of this paper is to solve generic 2-loop renormalization group (RG) equators 
exactly using the Lambert W function. Il| The Lambert W function has been introducec 

nn r 

previously to solve analytically the 2- and 3-loop RG equations for QCD.|2H4J (See also [5] 



for a review of various applications of the Lambert W function to QCD.) We apply the same 
function to solve generic 2-loop RG equations for theories with a Gaussian fixed point such 
as the O(N) non-linear sigma model in four dimensions. (This has been partially done in 

Bo 

In the following we wish to justify our purpose by reminding the reader of the generality 
of 2-loop RG equations. We consider two examples in four dimensions: QCD and the 4 
theory. Let us consider QCD first. 

Let A be the ultraviolet cutoff, and g\ be the bare gauge coupling, normalized appropri- 
ately. To construct the continuum limit (Ao — > oo) we must give a particular Ao dependence 
to g%: 

91 = lnAa + c l n l n Aa_l n Ak!) W 

where c = 6 ^ for QCD with no quarks. We can introduce a gauge coupling g 2 , renormal- 
ized at a renormalization scale /i, through the Ao independent constant in the denominator. 
If we choose 

then g 2 satisfies the 2-loop RG equation 

- ^/ = (g 2 ) 2 + c(g 2 r (3) 

exactly. 

We next consider the <fi theory defined by the bare action 

S = J d*x {^-d^ + ^ + ^ (4) 

with an ultraviolet cutoff Ao- This is a theory with the Gaussian fixed point = Ao = 0. 
We cannot take Ao all the way to infinity, but for Ao large compared with the physical mass, 
we obtain an almost continuum limit. For a given Ao, let the critical squared mass be 

m c 2 r (A ) = A A (X )A 2 (5) 



Then, to get an almost continuum limit, we tune the bare squared mass as 

(4tt) 2 



a 



C 



m 2 = A 4 (A )Ag + z m (\ ) [ ^ ] m 2 (6) 



3A 



where 
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Here, the renormalized coupling A is defined so that 



c = — , a = — - (7) 



The An dependence of z m and to is determined so that the theory with given m 2 and A have 
no An dependence except for non- universal contributions suppressed by inverse powers of A . 
The parameters A, m 2 renormalized at the scale fx satisfy the 2- loop RG equations (1-loop 
for m 2 ) exactly: 



- /7 A_3A_ _ _(_3X_Y _|_ C (-2L-Y 

-fi-^m 2 = aj^m 2 



(9) 



We have thus reminded the reader that renormalization schemes exist so that 2-loop RG 
equations become exact. Hence, solving 2-loop RG equations amounts to solving general 
RG equations. This paper is organized as follows. In sect. [Ill we solve the generic 2-loop RG 
equations (Q exactly in terms of the Lambert W function. (This has actually been done 
already in sect. II.B.l of Jfjl.) Then, in sec. we give the main results of this paper by 
constructing two physical parameters: one corresponding to the dimensionless coupling and 
the other corresponding to a physical squared mass. In sect. IIVI we invert the construction 
and express the renormalized parameters in terms of the physical parameters. In sect. |V] we 
generalize the exact effective potential in the large N limit of the O(N) linear sigma model 
7| to construct a trial effective potential for finite N, fully consistent with the 2-loop RG 
equations. 

In this paper we adopt the convention to fix the renormalization scale at fi — 1. Hence, 
a squared mass parameter acquires the canonical dimension 2 in addition to the anomalous 
dimension in its RG equation. 



II. GENERIC 2-LOOP RG EQUATIONS 



We consider the following generic 2-loop RG equation: 

d 



and 



—x = —x 2 + cx 3 
at 



-rrn 2 = (2 + ax)m 2 

(JjL 



For example, in the O(N) linear sigma model in four dimensions, we find 

9N + 42 



° (A + 8) 2 

for the self-coupling, and 

N + 2 
a = -NT~8 

for the squared mass parameter. 

In the following we assume c > and < x <C -. Let us define a mass scale by 



A(«) S {e±->(l-l)f>l 



This satisfies 

A(x) is of the same order as the UV cutoff. We can invert the definition of A(x) to 
x in terms of A(x). Since 

A(,)».,i-'(l-l) 

we obtain 

l-l = ^(A(x)^) 
where W is the upper branch of the Lambert W function defined by 

W(x)e w(x) = x 

for x > — |. (See Appendix 1.) 

We now define the running parameter x(t] x) by 



or equivalently by 

x(t; x) = — T r-- (20) 

c(l + w((e*A(x))-)) 

so that it satisfies both 

d t x(t; x) = —x(t; x) 2 + cx(t; x) 3 (21) 

and the initial condition 

x(0;x)=x (22) 



(Eq. (120]) agrees with (26) of [6] which gives the same result for c = 1.) 
Analogously, the running parameter m 2 (t; x, m 2 ) defined by 

i 



ri*~{t; x, m 2 ) = e 2t m 2 [ cS( f x) _ - ] (23) 



1 



c.r 



satisfies 

d t m 2 (t; x, m 2 ) = (2 + ax(t; x)) m 2 (t; x, m 2 ) (24) 

and the initial condition 

m 2 (0]x,m 2 ) = m 2 (25) 

III. PHYSICAL PARAMETERS 

We now apply the results of the previous section to construct physical parameters. We 
first observe that the combination 

2 -2c( J — l) 

is an RG invariant. We then note that any RG invariant of x and m 2 can be obtained as a 
function of the above RG invariant. 

Let x p h(x,m 2 ) be the RG invariant satisfying 

x ph (x, 1) = x (27) 

To obtain an explicit expression for Xp^x, m 2 ), let us write it in the form 

( x ,m 2 ) = f ( ( m 2 )-^J^- (—-l) e^r^- 1 )) (28) 
V a + 2c \cx J J 



x 



ph I 



where we assume m 2 > 0. The condition (1271) implies 

/ ( — -4r- f— — l) e ^(^ _1 ) I = ./■ 



a + 2c Vex 



(29) 



We can rewrite this as 



2c 



1 



(30) 



a + 2cyc/(se s ) J a + 2c\cx 

For small x < 1, we find s ^> 1 for a + 2c > 0, and — s ^> 1 for a + 2c < 0. The above 
equation is solved by the Lambert W function as 



2c 



1 



W(se s ) if a + 2c> 
iy_i(se s ) if a + 2c<0 



a + 2c \cf(se s 
We thus obtain the physical coupling 

a + 2c Vcx ph (x,m 2 ) / | ^ ^ m 2)-^_^_ (_L _ f) e » 



(31) 



if a + 2c> 



if a + 2c < 
(32) 

We plot the left-hand side as a function of m 2 for x = 0.1 assuming c = |= and a = — |. 



34/25 (27/(17 xph)-1) 



1 2 3 4 5 
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1. 



x p h increases 



FIG. 1. Plots of x ph (left) and §§ - l) (right) for x = 0.1, c = ffr,a - ... 
monotonically function of m 2 . It vanishes clS — j 2" 777^ —>• 0, and approaches -rX as to 2 — > oo. 



Though it is not obvious, the physical coupling x p h admits an asymptotic expansion in 
powers of x. This is because be defined by the differential equation 

d 



dt 



x 



P h 



(— x 2 + cx 3 ) Or + (2 + ax)m 2 d m 2 



X 



ph 







(33) 



and the initial condition ([27]) . For small x <C 1 , we can expand x p h asymptotically in powers 
of x in the form 



x p h(x, m 2 ) = x 



1 + x n p n (lnm 2 ) 



n=l 



(34) 



where p n is a polynomial of degree n satisfying p n (0) = 0. In principle, this can be shown 
directly from fl32|) . but it is more easily shown from the differential equation fl33|) . 
The physical coupling can also be given in the form of a running parameter: 



where t p h(a;,m 2 ) satisfies 



x p h(x, m 2 ) = x(— t p h(x, m 2 ); x) 



—t ph (x,m ) = 1 



and 



t ph (x, 1) = 



(35) 

(36) 
(37) 



To find t p h(x,m 2 ), we use the defining equality W(xe x ) = x to obtain 



i = w 



v CXph ( 



X Q CX P h( X > m ) 



VP A(z) 



_ i e c > (l ' m 1 -A(x)" 



VP(V^ h{a; ' m2) A(a;)^) 



Hence, we obtain 



t p h(x, m 2 ) = In A(x) — cln 



- 1 e cx p h 



cx h (x,m 



K ycx p h(x, m 2 ) , j 
In addition to the physical coupling, we can introduce a physical squared mass by 



This satisfies 



and the initial condition 



m 2 h (x, m 2 ) = in 



d 



2 / cx ph (x,m 2 ) 



- 1 



- 1 



—m ph (x,m ) = 2m ph (x,m 



™> 2 ph {x, 1) = 1 



(38) 
(39) 

(40) 

(41) 
(42) 



Using (1321 . we can rewrite the physical squared mass as 



W should be replaced by W-i if a + 2c < 0. The physical squared mass also admits an 
asymptotic expansion in x just as x p h- 



m h {x, m ) = m 



m 



n=l 



where q n is a polynomial of degree n satisfying q n (0) = 0. 




(44) 



FIG. 2. Plot of m 2 h and m 2 : We plot m 2 h for x = 0.1, c = i=, a 



P h 
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3" 



p h is essentially 

a monotonically increasing function of m 2 , even though it eventually starts decreasing when m 2 
reaches the UV cutoff scale. 



IV. x, m 2 IN TERMS OF x ph , mj h 

In the above we have introduced two physical parameters x p h, m, 2 ^ as functions of x, m 2 . 
We can invert their relations to express x, m 2 in terms of x p h, ^ph- We first rewrite (1321) and 
flU as 




where we assume m 2 > 0, and should be replaced by W-\ if a + 2c < 0. Substituting the 
second equation into the first to eliminate m 2 , we obtain 

(l_l).±-i = (m yi(_L _!).=*-' (47) 
8 



This gives 



1 „./, 1 



— -\ = W[ (ml h )-c ( — 1 ) e ca; ph " I (48) 



which is valid irrespective of the sign of a + 2c. Hence, we obtain 



11 

x = ? N (49) 

C l + W r ((m5 h )i(^-l)e-*- 1 J 



Using this result, we then obtain 



2 2 I ca; ph 



1 -1 



pM J_ _ i 



-p h { \ i a !> (so) 



V. A TRIAL EFFECTIVE POTENTIAL CONSISTENT WITH RG 

In jzj], the effective potential for the large AT limit of the O(N) linear sigma model in four 
dimensions has been obtained as 



2 



^Kfr(^) = mp h U m 2 + (4^) 2 x- ) (51) 



where v is the VEV of the scalar field with no anomalous dimension: 

Tt v = « < 52) 

In the large N limit, we obtain 

c = 0, a = -l (53) 



(54) 



so that 

A(x) = e* 
m 2 h (x,m 2 ) = A(x) 2 exp[H/_ 1 (-2^); 

In the symmetric phase m 2 > 0, m 2 h (x,m 2 ) gives the physical squared mass of the scalar 
fields (J) 1 (I = 1, • • • , N). In the broken phase m 2 < 0, the physical squared mass vanishes at 

-.2 



2 ~ m 

v = wri (55) 
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To generalize (15T|) for a finite N, for which c and a are given by (|T2l) and (fl3|) . we may 

try 

JrKfr(^) = ^(Xp h )m^ h ^ar, m 2 + #(x ph ) ^ - cj y^J > (56) 

which is a monotonically increasing function of v 2 . Here z and g are positive functions of 
the RG invariant 

x ' P h = x ph(x, \m 2 \) (57) 

which is well-defined irrespective of the sign of m 2 . Note that the term added to m 2 satisfies 
the same RG equation as m 2 : 



d 
dt 



1 



— c 



x 



a v 2 



(2 + ax) 



1 



— c 



X 



a y 2 



(58) 



There is no justification for (156]) except that it is fully consistent with RG, and that it gives 
the correct result in the large N limit where z and g are mere constants: 



z = 1, g= (4tt) 



(59) 



Note that in the broken phase m 2 < 0, the effective potential (or equivalently f l5"6"|) ) is defined 
only for 



2 \ 2 

V >v min 



-2m 2 



c) > 



(60) 

g{x ph ) \x / 

Since the right-hand side of (156]) is monotonically increasing with v 2 , the effective potential 
is minimized at v 2 = v^ in . 

The main advantage of the assumption ( 156]) is its integrability. To integrate ( 156]) with 
respect to v, we use (143]) to write (156]) as 
d 



— 2-KffO) = ^Oph) 
a 2 



m 2 , 



— c 



X < W 



2c 

2c a + 2c 



a + 2c \A(x) 



2 2 

m ( i \ v 



where W should be W_i for a + 2c < 0. Denoting 



i ■ 

' a + 2c 



(61) 



we can rewrite the differential equation for V e s as 



(62) 



&(^ph 



V 



' a + 2c 



(63) 
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We thus obtain 



2c 



^ z M Hx) ^ d H^ w {^^ ]} (64) 

Using the formulas 

POD 

/ dss^- 1 W{s) a 
= -{-p)- 1 -"-? [/3r (a + (3, -PW(s)) - T (a + /3 + 1, 

(,3<0,s>0) (65) 
J^dsi-sy-X-W-tis)) 

= -/3~ 1 -"-' 3 [/?r (a + /3, -0W_i(s)) - T (a + + 1, -/9W_i(s))] 
(/3>0,s<0) (66) 

where 

/■oo 

r(a,z) = / dtt a - l e- 1 (67) 
is the incomplete gamma function, we finally obtain 

Kfr ( v ) = 4442 2a+8c - 1 c 4c A(x) 4 

x {2(a + 2c)r (-a - 4c, 2(a + 2c)VU(s)) + Y (—a — 4c + 1, 2(a + 2c)W(s)) } (68) 

where 

2 2c 1 

S = C a + 2c 77 a + 2c (69) 

a + 2c ' K ' 

Note that W should be replaced by W_\ for a + 2c < 0. 



VI. CONCLUSIONS 

In this paper we have constructed two physical parameters x p ^ and by solving generic 
2-loop RG equations analytically in terms of the Lambert W function. In addition we 
have constructed explicitly a trial effective action, which is fully consistent with RG, by 
generalizing the analytic expression for the large N limit of the O(N) linear sigma model in 
four dimensions. 7|] The trial effective potential is, however, at best a wild guess at the true 
effective potential. Its only merit may be that it gives an intriguing example of what RG 
improved perturbation theory can produce. 

11 



The closed- form analytic expressions for x p h (given by (132]) ) and (given by ( 140]) ) sum 
the corresponding perturbative series. Further studies may elucidate the precise asymptotic 
nature of the perturbative expansions, as has been done for QCD. 2], 3]] 



Appendix A: The Lambert W function 



The Lambert W function is defined implicitly by 



W(x)e 



W(x) 



X 



(Al) 



or equivalently by 



W(xe x ) = x 



(A2) 



Restricted to real values, the function has two branches: the upper Wq(x) > —1 defined for 
x G [— e -1 , +00) and the lower W-\{x) < —1 for x G [— e _1 ,0). (See Fig. [31 ) For simplicity, 
we denote W as W in this paper. We obtain the following asymptotic expansions: 




FIG. 3. The real valued Lambert W function has two branches: the upper Wq (solid) and lower 
W-i (dashed). 



1. For x ^> 1, 



2. For -x < 1, 



Wn(x) = In x — In In x + O 



'minx" 
In x 



W-i(x) = ln(-x) - In (- ln(-x)) + O 



In (— ln(— x)) 
v ln(-x) _ 



(A3) 



(A4) 
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Appendix B: Asymptotic free theories 



Let us quickly summarize the applications of the Lambert W function to asymptotic free 



theories. 



3( For asymptotic free theories, the generic 2-loop RG equation is 

4-sc = x 2 + cx 3 (Bl) 
at 



For example, in QCD with n/ flavors, we find 

_ 6(153- 19n f ) 
° ~ 33 - 2n f 

and in the O(N) non-linear sigma model in two dimensions, we find 



N-2 

The scale parameter in this case is defined by 



(B2) 



C = (B3) 



A(x) = (e"^- 1 ( — + 1) ) (B4) 



cx 

Inverting this, we obtain 

x = p ) rrv (B5) 

-c{l + W(-(A(x))°)} 

Hence, the running parameter is given by 

x(t;x) = y~ rry (B6) 

-c{l + W(-{e*A(x))°)} 

which satisfies 

d t x(t; x) = x(t; x) 2 + cx(t; x) 3 (B7) 

and x(0] x) = x. 

The large t behavior of x(t; x) depends on the sign of c: 

• For c > 0, x(t; x) diverges as t — > t ma x, where t max is given by 

e" w =e c A(x) (B8) 

• For c < 0, x(t; x) approaches — as t — > +oo. 
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